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Abstract 

By reviewing the application of the renormahzation group to dif- 
ferent theoretical problems, we emphasize the role played by the gen- 
eral symmetry properties in identifying the relevant running variables 
describing the behavior of a given physical system. In particular, we 
show how the constraints due to the Ward identities, which imple- 
ment the conservation laws associated with the various symmetries, 
help to minimize the number of independent running variables. This 
use of the Ward identities is examined both in the case of a stable 
phase and of a critical phenomenon. In the first case we consider 
the problems of interacting fermions and bosons. In one dimension 
general and specific Ward identities are sufficient to show the non- 
Fermi-liquid character of the interacting fermion system, and also 
allow to describe the crossover to a Fermi liquid above one dimen- 
sion. This crossover is examined both in the absence and presence 
of singular interaction. On the other hand, in the case of interacting 
bosons in the superfluid phase, the implementation of the Ward iden- 
tities provides the asymptotically exact description of the acoustic 
low-energy excitation spectrum, and clarifies the subtle mechanism 
of how this is realized below and above three dimensions. As a critical 
phenomenon, we discuss the disorder-driven metal-insulator transi- 
tion in a disordered interacting Fermi system. In this case, through 
the use of Ward identities, one is able to associate all the disorder 
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effects to rcnormalizations of the Landau parameters. As a conse- 
quence, the occurrence of a metal-insulator transition is described as 
a critical breakdown of a Fermi liquid. 
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1 Introduction 

There exist several physical problems where the use of simple perturbative 
methods seems to be deemed to failure from the start. In these cases, 
already in the leading perturbative correction, one faces singular terms 
which render the perturbative expansion meaningless, no matter how small 
the initial expansion parameter is. The great success of the renormalization 
group (RG) has been to turn, at least in a number of highly non-trivial 
cases, this apparent failure into a powerful tool. The RG conceptual scheme 
^121 teaches us how singular terms, arising in the perturbative expansion, 
give rise to the power-law and scaling behaviors characteristic of the critical 
phenomena |HiEIHl!E!- This is the most successful use of the RG, which we 
only briefly recall in Sec. (21 In particular, we introduce the basic elements 
of the RG approach such as universality, scaling, relevant and irrelevant 
variables, and discuss the Wilsonj^l and the field-theoretic RG|3 Ej |S1, 
which can be seen as two different implementations of universality. 

The origin of a singular perturbation theory is understood within the 
(f"^ model for critical phenomena, where <p(r) is the field whose thermal 
average specifies the order parameter cpo !) . Indeed, the dimensionless ef- 
fective coupling constant which appears in perturbation theory in d spatial 
dimensions is u/\t\'^^^, where e = 4— d, u is the coupling constant describing 
the self-interaction of the fluctuations of the field ip, and t = T — Tc mea- 
sures the deviation of the temperature T from the critical temperature Tc. 
This leads to a divergence when approaching criticality {t ^ 0), for e > 
(d < 4) [for the corresponding power counting, leading to the expression 
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of the dimensionless coupling constant given above, see Fig. da)]. The 
RG sums the leading singularities into power laws for the order parameter 
(fio ~ \t\^ , for the susceptibility x ^ j which measures the linear re- 
sponse to its conjugate external field h, for the specific heat C ~ etc. 
Here, /3,7,a... are called the critical exponents, or critical indices. 

We will also turn our attention to different uses of the RG approach. 
One class of such applications is met when the critical behavior is uncon- 
ventional This may occur, for instance, when it is difficult to identify 
the order parameter (which can be a complicated object, or a composite 
operator) and the related broken symmetry. The problem of the metal- 
insulator transition (MIT) in a disordered interacting electron system is an 
important example of this class |10[ IllL I12L I13L 1141 1151 [TB] . In this case the 
symmetries of the physical problem may provide the guiding framework 
to identify the effective theory. These symmetries can be translated into 
Ward identities which establish relations among the various terms of the 
skeleton structure of the problem, simplifying the RG treatment. 

The use of the RG, however, is not limited to the critical behavior. It 
turns out to be extremely useful also in the stable phases. Indeed, there 
are cases in which a divergent perturbation theory is the result of relevant 
low- lying excitations in low-dimensional systems, despite the fact that the 
system is in a stable liquid phase of the matter. This is the case, e.g., of 
the Luttinger liquid in one-dimensional interacting Fermi systems |17| [see 
Fig. nib)], or of Bose liquids in the superfluid phase for d < 3 
[see Fig. ^c),(d)]. The very condition of phase stability imphes that 
the susceptibilities (e.g., specific heat, spin susceptibility, compressibility, 
etc.) are finite, and the exact cancellation of the singular contributions is 
controlled by Ward identities. Nonetheless, as we shall show, a track of the 
singularities is usually left in the power-law behavior of some one-particle 
physical quantities. 

In the following, we will touch upon all these examples, often stressing 
the key role played in practice by the Ward identities which implement the 
fundamental symmetries of the physical system at hand, and for the last 
two cases allow for a complete solution of the low-lying excitation problem. 
The rest of the paper is organized as follows. In Sec. |U we recall a few 
basic concepts about the RG and how Ward identities may be used to carry 
out the renormalization program. In Sec. |2| we discuss the application of 
RG to stable phases. First, in Sec. 13.11 we examine the case of interacting 
fermions in one [201121] and low dimensions (22 EHl El ESI • Then, in Sec. 
13.21 we analyze the problem of interacting bosons |26[ I27| . In Sec. 01 we 
focus our attention on the problem of disordered interacting Fermi systems 

|2HlEniEniEIlE21ESlElESllSniE71EHllSni 
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2 The renormalization group and the Ward 
identities 



2.1 General remarks on the group transformation 

The first aspect that must be appreciated by looking at the historical de- 
velopment of the RG is the change of perspective with respect to the re- 
duction scheme which, before its introduction, was at work in formulating 
a condensed-matter theory. The fifties and the sixties of the past century 
had seen the flourishing of the quasiparticle concept. According to this, 
an interacting system may be effectively described in terms of "almost" 
independent elementary excitations (quasiparticles). Implicit in this atti- 
tude is the privileged role reserved to the "approximate" solution of the 
dynamical problem and to the development of the many-body theory, with 
respect to the statistical aspects. In the late sixties and early seventies, the 
realization of the importance of the large-scale fluctuations in approaching 
the critical point shifted the attention onto the statistical problem |4l)[ I41| . 
The divergence of the correlation length ^ ~ 1^1*^, when t — + 0, implies the 
strong correlation of a large number of degrees of freedom within a region of 
an increasing linear size ^, and then invalidates the previous quasiparticle- 
based reduction schemes. The collective critical phenomenon does not arise 
as a simple superposition of single microscopic events. Through the action 
of ordering forces, the probability distribution of the collective variables no 
longer obeys the central limit theorem, according to which the correct nor- 
malization for an extensive variable is provided by the square root of the 
number of degrees of freedom. This violation, via an anomalous normaliza- 
tion of block variables, is the mathematical manifestation of the existence 
of correlations at all length scales, so that the crucial hypothesis of statis- 
tical independence becomes progressively not applicable getting closer to 
the critical point. 

The very fact that ^ —i- oo when approaching criticality, while invali- 
dating the previous reduction schemes, provides the crucial key for the new 
reduction scheme, i.e., universality PHIET] . 

There are two equivalent forms of universality. First of all, the mi- 
croscopic details which specify the peculiarities of each individual system 
become gradually irrelevant when approaching criticality, and the inflnitely 
large number of degrees of freedom involved is well accounted for by a small 
set of relevant variables. Once the proper choice of the relevant variables 
is made (e.g., t and h ot ip), one can change the parameters Q which spec- 
ify the microscopic details (unless they assume a value which changes the 
symmetry of the problem, thus becoming relevant), without changing the 
critical behavior of the system. Said in other words, systems which differ 
in the irrelevant variables share the same critical behavior • This condi- 
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tion translates into the invariance of the (singular part of the) free energy 
T ^ or equivalently of its Legendre transform F = J d'^r h^p — JF, 

^(i,/i;{CJ) = ^(<',/i';{Ca), 

r(t,<^;{0}) = r(t',<^';{ca), 

i.e., a suitable rescaling of the relevant variables fully accounts for a change 
in the irrelevant variables. 

The second form of universality relies on the fact that the unit length 
scale is irrelevant when ^ ^ oo. This is translated into the statement that 
one can eliminate the degrees of freedom at short distance, by grouping 
microscopic variables into "block variables" , within an iterated procedure 
to build larger and larger blocks, without changing the critical behavior of 
the system |41| . 

The field-theoretic RG equations generalize the universality relations 
in the sense that they relate one model system to another by varying the 
coupling constant C and suitably rescaling the other variables and the cor- 
relation functions [?]■ In ordinary critical phenomena, described by the 
(p^ theory, ip (wave-function) and t (mass) must be rescaled. When the 
coupling reaches its fixed point, a two-parameter scaling follows. 

The Wilson new RG implements the second idea of universality illus- 
trated above, by eliminating the short-wavelength fluctuations of the field 

i.e., those with characteristic momenta A/s < g < A, where A is the 
ultraviolet cutoff, and s is the scaling parameter |2] . The real space version 
of the Wilson RG implementation is the bloc variable transformation!^ 

namiiini. 

The calculation of the critical indices was performed either numerically, 
through recurrence equations, or analytically, within the e expansion 14 7| . 
This latter procedure relies on the perturbative expansion of the RG trans- 
formation in the parameter e = 4 — d, starting from the known case d = 4. 
The direct calculation of the divergent critical quantities is thus avoided. 

In general, a sets of parameters {/ij} specifies a given Hamiltonian iJ 
in the space H, of the Hamiltonians. The action of the RG transformation 
T^s (which depends on the scaling parameter s) in this space preserves the 
functional F, so that symbolically 44 , 45 : 

F[7^.[i^]] = F[ff]. 

The above equation can be written in a differential form as 

dnV = 0. 

The functional derivatives of F yield the group equations for all the physical 
quantities. 
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A fixed point of the transformation is such that H* = TZs[H*]. At 
criticahty ^ = 00 and it remains infinite under iteration of RG transforma- 
tions. By iteration of the group transformation on the critical surface one 
model Hamiltonian is transformed into another, until a fixed-point simpli- 
fied Hamiltonian is reached as s ^ 00, which is manifestly scale invariant, 
and all the irrelevant transient terms have been eliminated. The relevant 
parameters define a finite set of relevant directions of escape from a given 
fixed point |3S1) and their scale dependence gives a microscopic definition 
of the critical indices. Universality corresponds to the domain of attraction 
of a given fixed point. Within a given domain, the transformations become 
asymptotically eouivalent [431 IH| . 

Standard renormalizability implies that the removal of the ultraviolet 
cutoff, A ~ 1/a — > 00 (where a is the lattice spacing or some other mi- 
croscopic characteristic length scale) only requires the definition of a finite 
(small) set of renormalized parameters. Relevant parameters and renor- 
malized parameters coincide asymptotically in the infrared region. 

Here a general problem arises. Recently, a renewed attention has been 
devoted to the exact RG equations 0^1 ■ In practice, however, one has 
to truncate these equations or reduce the number of flow parameters. In 
this last case, one has to initialize the action to follow a renormalized 
trajectory with a small number of flowing parameters, and an infinite set 
of parameters in the action must be specified. This specification is a highly 
non-perturbative condition. 

In practice, therefore, one should rather choose the best simple action, 
guided by the symmetries of the problem and by the physical conditions. 
The knowledge of the fundamental symmetries inherent to each specific 
problem, and in particular those of the order parameter, have to be assumed 
in order to make the proper choice of the basic variables on which the RG 
transformation acts. 

2.2 Various exploitations of the Ward identities 

The RG approach takes enormous advantage from the explicit implementa- 
tion of the symmetries of the original problem, which translate into Ward 
identities relating various independent quantities. 

First of all, the Ward identities allow for a reduction of the renormaliza- 
tion parameters. The case of quantum electrodynamics is a textbook exam- 
ple. The local C/(l) gauge invariance of the Lagrangian, which is invariant 
under the simultaneous transformation of the phonon field ^ + d^-d 
and of the fermion field ip ipe''^'^ ~ {l + ied)ip, where 1!) is the gauge field, 
allows to relate the self-energy and vertex renormalizations, whose correc- 
tions correspond to superficially divergent diagrams within perturbation 
theory in d = 3 + 1 dimensions. 
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Secondly, the Ward identities provide the guiding framework for the 
identification of the proper running variables (effective coupling constants), 
and how they are related to physical quantities. We will consider the case 
of disordered interacting electron systems. As we shall see, the various 
interaction amplitudes are dressed by disorder. The renormalization pa- 
rameters of the corresponding field-theoretical formulation of the problem 
(the non-linear tr-model |2HI^5) are identified in terms of physical quan- 
tities of the Fermi-liquid theory by requiring the gauge invariance of the 
skeleton structure of the response functions ^1 1^ |32 ESI IHH ■ 

Finally, the Ward identities control the exact cancellation of the sin- 
gularities in the various response functions within stable phases, despite a 
singular perturbation theory. In these cases, indeed, specific symmetries 
are related to additional Ward identities, which enforce the cancellations. 

The combination of the three procedures indicated above leads to the 
asymptotically exact description of: the Luttinger liquid in d = 1 [T71 1^ . 
and its crossover to the Fermi liquid in d > 1 ; the non- Fermi-liquid 

behavior of an interacting electron system with singular forward scattering 
in d > 1 1^ the Bose liquid in the presence of the Bose-Einstein 
condensate [5^; the disordered interacting electron systems |^ |^. In 
Tab. ^we provide a scheme of the symmetries and Ward identities which 
are relevant in the various physical problems discussed in this paper. 

Each specific system has been discussed at length in several papers. 
We emphasize here the general common aspects and give a unifying view 
of these apparently very different theoretical problems. 

3 Renormalization group and stable phases 

As we have already indicated in Sec. ^ there are cases where the pertur- 
bation theory is singular even in a stable phase. We first recall here the 
case of the Luttinger liquid in one dimension, its dimensional crossover to a 
Fermi liquid as soon as d > 1 in the presence of short range forces, and the 
non-Fermi-liquid behavior when the forces among the fermions are suffi- 
ciently singular. We then consider the low-lying excitations from a ground 
state of bosons in the presence of interaction, where the condensate changes 
the power counting and leads to singular terms in perturbation theory for 
1 < d < 3. The use of the RG, in these cases, acquires a different meaning 
with respect to critical phenomena, as the thermodynamic stability implies 
here a cancellation of singularities at any order in perturbation theory in all 
the susceptibilities, like specific heat, compressibility, etc. Behind these ex- 
act cancellations there must be specific symmetry properties, which can be 
used in the form of Ward identities to close the hierarchy of RG equations 
and solve the problem. Although a stable system is characterized by finite 
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response functions, the re-summation of singularities may still manifest it- 
self with an anomalous dimension as in the one-particle fcrmion Green's 
function. 

3.1 Interacting fermions in low dimensions 

Within a RG approach, the effective low-energy Hamiltonian which de- 
scribes a metallic system can be obtained via the Wilson-like iterative 
elimination of states outside a shell of thickness A around the Fermi sur- 
face lEI 1^31 ISU] . As A is iteratively scaled to zero a fixed point is reached 
in d = 3 which corresponds to a gas of quasiparticles with finite resid- 
ual Hartree-like interaction, described by the Landau functional. Ordinary 
three-dimensional metals are thus well described by the Landau Fermi- 
liquid theory |51| . The interaction is effectively taken into account in the 
various response function through a small set of Landau parameters, enter- 
ing the physical quantities, e.g., the specific-heat coefficient 7 = Cy/T, the 
spin susceptibility Xi the compressibility 9n/9/i, and the spectral weight 
of the Drude peak in the optical conductivity, all of them being finite. 

A finite wave-function renormalization Z < 1 reduces the discontinuity 
of the occupation number in momentum space at the Fermi surface with 
respect to the value of the Fermi gas {Z = 1). The elementary excitations of 
the system are the single-particle excitations (quasiparticles), which carry 
both charge and spin, and the collective charge and spin modes, which get 
overdamped when they enter the particle-hole continuum of excitations. 

In recent times, the interest for non-Fermi-liquid metallic phases has 
been rekindled by the observation of anomalous features in single-particle 
and transport properties of the cuprates, which have been interpreted as a 
signature of non- Fermi-liquid behavior '52'. These compounds are layered 
(quasi-two-dimensional) materials, which are insulating when stoichiomet- 
ric, and become metallic upon chemical doping. The metallic phase is 
strongly anomalous at low doping and gradually (and smoothly) evolves 
towards a Fermi-liquid-like metal at large doping. 

The breakdown of the Fermi liquid could be produced by the suppres- 
sion of quasiparticle spectral weight at the Fermi surface due to the gener- 
ation of an interaction-induced anomalous dimension in the wave-function 
renormalization, Z ^ Ip ~ Pf|'', where p — pp measures the deviation of 
the quasiparticle momentum p from the Fermi momentum pp. In such a 
case, Z vanishes as p — > pp, the low- lying single-particle excitations are 
completely suppressed, and the low-energy behavior of the system is domi- 
nated by the charge and spin collective modes. Since these have, in general, 
different velocities, the system is characterized by the so-called charge and 
spin separation. 

This behavior is achieved in d = 1, where the metallic phase is a Lut- 
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tinger liquid. The hint for the breakdown of the Fermi-hquid picture is 
given by the appearance of logarithmic singularities in perturbation theory 
as p — > [see Fig. ^b)]. These singularities are controlled by the RG, 
and an effective field theory results, with marginal forward scattering, the 
so-called Luttinger model, which is exactly solvable. Additional conser- 
vation laws, which are peculiar to one-dimensional systems with forward 
scattering, constrain the model and lead to a closure of the equations of 
motion. 

Indeed, the total charge and spin conservation is translated into the 
usual Ward identity 



which connects the irreducible density (Aq) and current (A) vertices to the 
electronic Green's function G. In addition, in this case, charge and spin are 
separately conserved at each Fermi point p = ipF, leading to an additional 
Ward identity [201 stating the proportionality of the current vertex to the 
density vertex, due to the unique momentum direction. 



where is the Fermi velocity and the direction versor p here (d = 1) is ±1. 
By means of the two Ward identities ^ and jSJl one is able to express the 
density vertex Aq as a function of the fermion propagator G, thus closing 
the Dyson equation of Fig. [S] 

The effective interaction D (which is represented by the tick-dashed 
line in Fig. O contains the charge and spin collective modes via the RPA 
re-summation of the bare interaction with the fermion polarization bubble. 
The RPA re-summation is shown to be exact by the use of the same Ward 
identities and (0). As a result density and spin response functions are 
finite, leading to a stable metallic phase. 

On the contrary, from the exact evaluation of G one finds that its behav- 
ior is controlled by a non-universal anomalous dimension rj. Single particles 
hence move as composite objects due to the strong mixing with the charge 
and spin collective modes, induced by the effective interaction. 

The question then arises whether the Luttinger-liquid behavior can be 
extended to higher dimensions. The analysis of the dimensional crossover 
shows that this is not possible if the interactions are short-ranged. Indeed, 
the generic integrals in d dimensions 



where Sd-i is the solid angle in d — 1 dimensions, are strongly peaked at 
= 0, TT for 1 < c? < 2, so that all the relevant vectors are still parallel 



wAo-q-A = G ^(e + cj,p + q) - G ^(e, q). 



(1) 



A = vppAo 



(2) 
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or antiparallel. In this case, the additional Ward identity l(2Jl is still valid 
asymptotically near the Fermi surface. However, the mixing with the col- 
lective modes is reduced as soon as d > 1, since the effective interaction 
D{q,uj) is averaged over the d — I momentum components perpendicular 
to the Fermi momentum. The effective interaction, which is marginal in 
d = 1, scales to zero for d > 1, and the system is a Fermi liquid. 

However, if we consider a long-range bare interaction among the fermions, 
V{q) ~ the effective interaction dressed by the RPA series is 

and is dominated by a collective mode u){q) — cg^^"/^, which is propa- 
gating and gapless for a < 2 (the parameter c depends on the microscopic 
parameters of the fermion model, and does not acquire singular corrections, 
see below). This singular behavior as q — > compensates the rescaling to 
zero of the effective interaction due to its averaging over the transverse 
momentum and leads to a non- Fermi- liquid behavior for a > 2d — 2 |24| . 

A singular interaction may appear in a system close to an instability, due 
to the coupling of the fermion quasiparticles with the critical fluctuations. 
Various proposals for the breakdown of the Fermi-liquid picture in the 
cuprates rely on the existence of a quantum critical point in their phase 
diagram 53_. 

The problem of fermions interacting via a singular interaction is plagued 
by infrared divergences and requires a proper RG treatment, supported by 
the Ward identities |25j . 

The vertex Aq induces a coupling constant g which in principle con- 
tains the divergences of the vertex (which requires the renormalization 
Zv), of the fermion propagator (renormalized by Z) and of the collective- 
mode propagator D (renormalized by Zo)- The coupling is accordingly 

1 /2 

renormalized as go = gZZjj jZy- However, the total charge conservation 
translated into the usual Ward identity J^lj gives, in the dynamic limit, 

ZZ^ = lim lim Ao = d^G'^ = Z'^ . (4) 

Moreover, since small- q scattering dominates, the charge is asymptotically 
conserved at each point of the Fermi surface, leading to the additional Ward 
identity ||2J). When the vertex Aq is accordingly expressed in terms of the 
fermion propagator G, one realizes that the effective interaction Z?, i.e., the 
collective-mode propagator Q does not get anomalous corrections beyond 
the RPA approximation, i.e., Zd = 1. In fact, the vertex corrections 
to the bare polarization bubble asymptotically cancel with the selfenergy 
corrections, thank to the additional Ward identity ||2Jl. 
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Therefore, together with Eq. this implies that g — go, and the 

dimensionless running couphng constant u — g^js^^, where s is some in- 
frared cutoff which plays the role of the scaling parameter, evolves under 
RG with its bare dimension a;„ = l — (i + a/2ass^0. Therefore at the 
dimension d = (ic = l + a/2 the coupling is marginal (xu = 0). For d < dc 
the bare dimension Xu is positive, the theory scales to strong coupling and 
the Fermi liquid breaks down. Finally, for d > dc, Xu is negative and the 
effective coupling scales to zero. Thus Z stays finite and the system is a 
Fermi liquid. 

From the one-loop perturbative RG, one finds the explicit form of the 
wave- function renormalization as a function, e.g., of the quasiparticle en- 
ergy e. At d = dc, Z ^ |£|"°, vanishes with a non-universal exponent 
Mo ~ (?o, as in the Luttinger liquid. For d < dc, instead, Z scales to zero 
as a stretched exponential, for e ^ 0. In both cases, the single-particle 
spectral weight is suppressed as the Fermi surface is approached. 

3.2 Interacting bosons in low dimensions 

The theory of the interacting Bosc system has been motivated by the ob- 
servation of superfluidity in liquid helium for T < Tx — 2.17K. The connec- 
tion between the form of the spectrum of the low-lying excitations and the 
superfluid properties of the system was established by Landau via the cri- 
terion which requires a finite minimum slope of the quasiparticle dispersion 

£(q) m 

The simplest interacting Bose system assumes a quartic contact interac- 
tion V for the boson field (p. The mean-field solution due to Bogoljubov j55| . 
is obtained by linearizing the interaction term via the factorization of the 
condensate density nj^^ — (if) — (<^^)- The free-particle spectrum q^/2m 
is converted at low momenta into a linear spectrum e{q) = y^2novq = coq, 
characteristic of a sound-like excitation in the liquid, leading to superflu- 
idity. 

However, the Bogoljubov solution leads to a non-zero anomalous self- 
energy E — v{(p(p) — vTiQ. Moreover, the flrst corrections to the Bogoljubov 
approximate solution, within the so-called pairing approximation, lead to 
a spurious gap in the excitation spectrum |[56j. Both these results are 
in contrast with exact results, which show that no gap is present in the 
spectrum (Hugenoltz-Pines theorem) |S7], and that the anomalous self- 
energy vanishes, E = j^Hl- 

On the other hand, standard perturbation theory is plagued by infrared 
divergences in d < 3, due to the Goldstone sound mode |18l I19| . despite 
the fact that the superfluid phase is a stable liquid phase of matter. All 
these aspects call for a cautious analysis of the problem, as it was first rec- 
ognized in Ref. J^l- Benfatto used the Wilson RG approach to determine 
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the scaling behavior for a number of running variables |27]. However, his 
treatment was limited to d = 3 and did not take advantage of a system- 
atic implementation of the Ward identities. We summarize the results in 
generic d dimensions, fully exploiting the Ward identities, along the lines 
ofRef. 

Within the standard formulation of the problem, an external source 
Ay = {n{r), A(r)), — 0, ...,d, coupled to the boson density (current) and 
a field h{r) conjugate to (p{r) are introduced to generate the correlation 
functions. The original problem is recovered when A — > 0, ft. ^ and 
II, taken as a constant, is identified with the chemical potential of the 
liquid. By introducing the longitudinal and transverse components to the 
direction along which the symmetry is broken, we have tp — ipL +^^t, with 
(vl)?!— ►o = y/^i a-iid h — Hl + ikx- Functional derivatives of the free 
energy T[A^,hi] with respect to the external sources produce the various 
density, current expectation values and Green's functions, e.g., ipio{r) = 
{(fi) = 5J-/5hi, with i = L,T. 

The Legendre transform r[Ay , t^io] = / d'^T^ X^i ^i'/'jo — ^[Av, hi] is the 
generating functional of the various vertex functions, e.g., hi{r) = 5T/5ipiQ. 

In this representation the mean-field Green's functions in frequency and 
momentum space read 



Gll 
Glt 



9 9 ' 

Cor 



2 9 ' 



(5) 



G 



TT 



9 9 ' 



whence it is seen that the most singular contribution in the infrared is 
carried by Gtt- In the presence of the condensate, the four-leg vertex 
associated to the coupling constant v yields various four- and three-leg 
vertices, associated to the various coupling constants vllll, etc., vlll, 
etc. (see Fig. Starting from Eq. the dimensional analysis, in 

units such that [q] ^ [u] — 1, identifies nine relevant or marginal running 
variables with singular corrections in d < 3. 

However, the number of independent parameters is reduced by the ex- 
ploitation of the local gauge invariance of the model. The free energy JF, 
or equivalently F, are invariant under the local transformation generated 
by the operator 



r 



cost? — sini? 
sin d cos "d 



1 + d 



which acts in the L, T space, i.e.. 



(6) 



13 



with i — L,T and j — L, T. 

Functional derivatives of T with respect to i?, A^, and (p generate the 
various Ward identities. The first two are obtained by a functional deriva- 
tivation with respect to 'd and then either with respect to 1^9^ or ipx- In 
this way the Ward identities connect the two-point vertices to the external 
fields, showing, e.g., that when h ^ no gap is present in the spectrum. In 
this way the Hugenholtz-Pines theorem is simply recovered. Three-point 
vertices are related to two- and four-point vertices via three relations, which 
are the implementation of the continuity equation for this special case in 
the presence of the condensate. By then, all the strongly relevant coupling 
constants are fixed, and one is left with four marginal running variables. 
However, three of them are connected to physical quantities. In fact, they 
are related to the derivatives of the vertices with respect to frequency or 
momentum, and the corresponding Ward identities involve density and 
current correlation functions, thus implying relations to superfluid density, 
condensate compressibility, sound velocity, all of them free of divergences. 
Therefore, these identities guarantee the exact cancellation of the singular 
contributions in these variables, which are then RG invariants. 

After all the identifications are made, one running variable is left, e.g., 
the longitudinal two-point vertex coupling v^l. 

The one-loop RG equations reproduce this situation of being left with 
only vll as running coupling, with vll ~ s'^, with e = 3 — d > 0, at d < 3. 
In this case the anomalous self-energy E vanishes. 

For d > 3, instead, vll ~ '2fJ- is finite, and the Bogoljubov result is 
recovered. In both cases, the spectrum is linear, leading to sound-like 
low-energy excitations. However this spectrum is realized in a completely 
different way in d < 3. 

As in the Luttinger liquid for the one-particle Green's function, the 
behavior vll ~ s'^ is fixed exactly by the Ward identities, which allow to 
close the (dominant part of the) equation for vll (see Fig. since the 
three-leg vertex vttl is identified with vll/vlo via Ward identities. When 
dealing with RG loop expansions, further corrections to the one-loop RG 
results cannot change the power-law exponents. 

4 An example of unconventional criticality: 
the metal- insulator transition 

4.1 Metal-insulator transition in disordered electron 
systems 

When we deal with disordered electron systems, the singularities do not 
come directly from integrations of loops in terms of the original fermion 
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propagators, as in Fig. ^ but rather from the soft collective modes related 
to diffusion, in terms of which an effective action can be derived for both 
non- interacting [591 l6(Jj and interacting systems [2H|- In this case, we will 
therefore enter in rather more details with respect to the previous topics. 

According to basic quantum mechanics, within the one-electron approx- 
imation, insulating and metallic behavior occur as a result of the complete 
or partial filling of the highest occupied energy band of a solid, respectively. 
Beyond the one-electron approximation, correlation effects may lead to a 
drastic rearrangement of the energy levels and give rise to insulating be- 
havior in situations where the metallic one is expected on the basis of band 
filling. This is usually referred to as a Mott phenomenon. 

Scattering from randomly located impurities may also lead to a transi- 
tion from a metallic to an insulating behavior. In this second case, rather 
than a rearrangement of the energy levels, the coherent scattering from 
the impurities modifies the phase of the electron wave function in such a 
way that it may result localized in space, thereby yielding an insulating 
behavior |61) . Experiments were performed both in metallic and semicon- 
ducting systems. These latter, in particular, turned out to be better suited 
for the study of the MIT, since the degree of disorder could be changed by 
the doping level, as it happens, for instance, in silicon doped with phos- 
phorus (Si:P). There exist several review articles, which summarize the 
development of the field at various stages ^1 ^fl E", 13 , 1^ El ■ 1^ 
recent years, there has been a renewed interest in the MIT, which is very 
actively investigated in the two-dimensional electron gas of MOSFET and 
heterostructure devices (for a recent review see, e.g., Ref. 

At theoretical level, a crucial element was the discovery that the quan- 
tum corrections to the classical Drude formula for the electrical conductiv- 
ity give rise to singular terms in low-dimensional systems. These corrections 
arise as a result of the quantum interference of electron waves in disordered 
systems. There are two types of such terms. The first, known as the weak- 
localization (WL) correction, is a purely one-particle effect and is due to 
the interference of time reversed trajectories [SSI IM] (see Fig. 0J). The 
second is a consequence of the enhancement of the electron-electron inter- 
action in a disordered system and is usually referred to as electron-electron 
interaction (EEI) correction [^31201 (see Fig. EJ. 

In two dimensions, in particular, both the WL and the EEI correction 
to the conductivity are logarithmic, 

A<TWL = j-\n{L^/l), (7) 

TT n 

AaEEi = ^HLt/1). (8) 
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Both the above effects have been discussed in great detail in the hterature 
and we do not provide a derivation of Eqs. (|7I8|I here, but hmit ourselves 
to a comment on the EEI correction, after Eq. H12|l . 

In Eqs. H7I8|I . the argument of the logarithmic corrections contains the 
ratio of two length scales. The first is the mean free path, I = vpT, r being 
the elastic scattering time, which sets the microscopic scale, beyond which 
the system behaves diffusively. The second length scale differs in the two 
cases. For WL it is L^, the scale over which inelastic scattering starts to 
destroy the interference effects and one can show that it becomes infinite at 
zero temperature |67| . For EEI it is usually given by the thermal length Lt ■ 
In a diffusive system, all length scales correspond to characteristic times via 
the diffusion constant D = VpT/2: the elastic scattering time t — P/D, the 
dephasing time — L'^/D, and the thermal time tt = fi/kBT — L^/D. 
As a result, the above corrections become singular at low temperature. 

By recalling the Drude formula for the electrical conductivity a — 
2e^NoD, where A^o is the single-particle density of states at the Fermi 
level, one sees that the corrections H7I8|I in units of a are controlled by the 
factor t — \ / {AiT^ NqDTi) (not to be confused with the deviation from the 
critical temperature of the previous sections). The parameter t is related 
to the two-dimensional conductance g = G'/(e^/ft.), by t = l/(g7r). In good 
metals, normally, g S> 1 and t is a small parameter. 

The WL correction has been used to formulate a scaling theory of the 
MIT [531221 (fo'^ * pedagogical introduction see, e.g., Nagaoka's and Kawa- 
bata's contributions in Ref. ^Jj)- This theory has a unique scaling vari- 
able, the dimensionless conductance g, which according to the Thouless's 
argument describes how the eigenstates of the system change when we add 
together blocks of size L to make blocks of size 2L |^. There are two en- 
ergy scales to be considered. The first is the energy level spacing that goes 
as Ai5 = l/{dN/dE) = 1/{NqL'^). The second is the energy perturbation 
due to the change in the boundary conditions. This may be related to the 
travel time in a diffusive system of size L, i.e., 5E Tl/t — Ti^J D jL?- . The 
ratio of these two energies gives 

— ^ ^ ^ Td-2 ^ 9_ 

AE 27r2e2//i 2tt' 

Starting from the metallic side, where tr is a constant, we see that in two 
dimensions the conductance g becomes marginal. By assuming a scaling 
behavior, one can derive a RG flow equation for g from the WL correc- 
tion in Eq. Q. As a result one predicts that, in two dimensions, the RG 
flow drives the system to an insulating ground state, with g ^ 0. The 
EEI shows qualitatively the same effect and one wonders whether one has 
to introduce further scaling parameters to take into account the interac- 
tion strength 69 . The answer to such a question is far from trivial and 
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has required quite some work. In particular, it has been understood that, 
within a Fermi-hquid description, the Landau parameters provide the ad- 
ditional running variables whose RG flow has to be studied together with 
the conductance g |2H1IS- In order to emphasize the physical origin of the 
additional running variables, we show in the next section how the interplay 
of interaction and disorder gives rise to singular corrections to thermody- 
namic quantities. This is to be contrasted with the WL correction which 
leaves the thermodynamic properties unaffected. 

4.2 Thermodynamics 

We present the derivation of the correction to the thermodynamic poten- 
tial due to the interaction and disorder PUIEHIED]. The coupling to the 
impurities is described by the Hamiltonian 



Greek indices indicate spin indices and sum over repeated indices is under- 
stood. To begin with, let us consider the first-order exchange-interaction 
correction to the thermodynamic potential (see Fig. O 



where we have introduced a retarded electron-electron interaction ^(x — 
x',w) to take into account screening effects which we will come back to 
later on. lo and e are Matsubara boson and fermion frequencies, respec- 
tively. Notice that the electron Green's function depends separately on 
its spatial arguments, since translational invariance does not hold in the 
presence of disorder. In general, however, the impurities are randomly dis- 
tributed and it is enough to average over them. For the average quantities 
the translational invariance is restored. The average can be performed by 
means of the well-known impurity techniaue|54|. In Eq. (O this implies to 
average the product of the two Green's functions. We assume a Gaussian, 
5-correlated impurity potential U{y) with 




while the electron-electron interaction is of the standard form 





(9) 



{U{v)U{v')) 



1 



5(r-r'), 



(10) 



2ttNot 
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and use a self-consistent Born approximation for the electron self-energy 
(see Fig. 

E(r,i;r',i') = ^^^G{r,t;r,t')S{r - r'). (11) 

When the above self-energy is used to evaluate the Green's function, one 
obtains a self-consistency equation which leads to a finite lifetime r. The 
average procedure is limited to correlations of impurity insertions according 
to Eq. H10() . Insertions made on the same Green's function are automati- 
cally taken into account in Eq. (|ll|l . When one has an average of two or 
more Green's functions, as in evaluating the response functions, correlations 
between insertions on different Green's functions must also be considered. 
To this end one re-sums the infinite series of the so-called ladder diagrams 
(see Fig. |7|l 

i(q,w) = 1 + + + . 



1 



where 



^=^EG(P,^)G(p + q,e + C.) 
P 



The sum over the momentum p differs from zero only when the Matsubara 
energies e and e + lo have opposite signs. For vanishing frequency and 
momentum, — > 1 signaling the emergence of a (diffusive) pole in the 
ladder. Indeed, in the low-frequency and momentum approximation, lut <ti 
1 and vpqT <ti 1, which defines the diffusive regime, the ladder sum reads 

L(q,uj) = ^—^ , , . (12) 

The above equation is telling us that the two-particle propagation has a dif- 
fusive form. This means that, whereas the single-particle Green's function 
is exponentially decaying over a mean free path, I = vpT {vp is the Fermi 
velocity), the collective density fluctuations propagate diffusively over large 
distances. It is the long-range diffusive character of the density fluctuations 
which is responsible of the logarithmic behavior of the EEI corrections to 
the conductivity and to the thermodynamic properties in two dimensions, 
as we are about to see. Technically, the logarithmic corrections arise from 
the integration over the diffusive pole. Each momentum integration, in two 
dimensions, yields a factor 1/D ~ t, the expansion parameter. We would 
like to emphasize that the emergence of critical massless modes is highly 
non-trivial and shows how the construction of the effective action is in this 
case rather unconventional, to say the least. 

The diffusive form of the two-particle propagator (usually called diffu- 
son) makes the small-g region more relevant. As a consequence, the Fourier 
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transform of the interaction in Eq. becomes important only for small 
momentum transfer, i.e., V{q = 0). On the other hand, the corresponding 
Hartree diagram contribution, after the average over the impurities, selects 
the interaction at large momentum transfer, i.e., V{q = 2pf). 

The good-metal condition, g ^ 1, can be written in an equivalent form 
as Epr/h » 1, which implies that the disorder only affects states within a 
small shell h/r away from the Fermi surface. Under these circumstances, 
since disorder modifies only states near the Fermi surface, interaction effects 
at larger energy can be taken into account via the Fermi-liquid theory, 
which goes beyond the first-order perturbation theory in the interaction, 
by replacing V{q = 0) and V{q = 2pf) with the corresponding Fermi- 
liquid scattering amplitudes Fi and whose lowest order diagrams are 
depicted in Fig. |S| Since, in the absence of spin-flip mechanisms, the 
total spin entering the two-particle propagator is a conserved quantity, it 
is convenient to introduce the singlet and triplet scattering amplitudes 

r. - Fi - Ifs, F, ^ ^T2. 

We recall that the scattering amplitudes Fs and Fj are related to the Lan- 
dau Fermi-liquid parameters F^^^ which enter in the compressibility and 
spin susceptibility, by 

1 FO 1 F^ 

F, = — Tt = —. 



2No 1 + FO' 2Nol + FO ' 

From now on, when necessary, A^o is assumed to include the Lan- 
dau effective-mass correction. The scattering amplitudes are dynamically 
screened by the diffusive density fluctuations, leading to dressed scattering 
amplitudes 

r ( ^-v + 1^1 

- ^-*Z?g2 + (lT2iVoF,,)M- 
However, this dressing does not alter the degree of divergence of a given 
diagram. 

We are now ready to present the expression for the correction to the 
thermodynamic potential. By means of the previous substitutions, this is 
obtained to all orders in the interaction, but to first order in the expansion 
parameter t. It is convenient to use the standard trick j71| of multiplying 
the interaction by a parameter A and integrating over it between and 1. 
We then get 

Dq^ + il-X2NoTs)\to\ 
NoTM 



An= - J dX 



E 



Dq^ + (1 + X2NoTt)\Lu\ - iM{l + 2NoTt)LUsSgn{Lu) 
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To allow for the calculation of the spin susceptibility, in the above equation 
we have introduced also a magnetic field via the Zeeman coupling, ujs = 
gfJLBB and M labels the triplet states. Notice the introduction of the 
Fermi-liquid renormalization of the Zeeman energy |72| . After evaluating 
the integrals one gets 



where 



Alio = -tNaT^{2NQTs - QN^Tt) \n{TT) + A 

Aril = tNoNoTt[l + 2NoTt)ul\n{uo,T), 
m^ix) = {\ + 2NoTt)I{x)-I{{l + 2NoTt)x)~AANoTt, 

poo 

I{x) = / dyb{y)[{v-x)\n\y~x\ + {v + x)\n\y + x\], (14) 

where h{x) is the Bose function and A « —0.24. Eqs. H13I14() are orga- 
nized to show the zero- and the large-magnetic-field behaviors, which are 
connected by the crossover function Ari2(a;). Notice that at small x, one 
has /^Vl2{x) « — A'ort(l -I- 2iVor()a;^ Inx. As a result, the corrections to the 
specific heat and to the spin susceptibility read 

5Cv = Cy,oi(^or.-3iVort)ln(TT), 
5x = -Xo4tA^ort(l + 2iVor,)ln(TT), (15) 

where Cvfi — 2tt'^NoT/?> and xo = ^^0(5/15)^/2 arc the non-interacting 
values. We note that there is no dependence on the chemical potential 
implying that there are no singular corrections to the static compressibility 
dn/dfji, i.e., 5{dn/d^) — 0. As a result, the compressibility has the Landau 
value 

1^ = TVo (1 - 27Vor.) ^ N^Zl (16) 

Eg. ljlSI) contains the leading divergent perturbative terms to thermo- 
dynamic properties of a disordered Fermi liquid. Before proceeding to the 
next sections, where we discuss the renormalized perturbation theory and 
how to obtain the equations for the RG flow including transport properties, 
we remark that in various physical systems the above corrections can actu- 
ally be observed already in the metallic phase, where they are still within 
the reach of perturbation theory jl3j . 

4.3 Ward identities and renormalized Fermi liquid 

In order to get the RG equations, one has to absorb the logarithmic cor- 
rections in terms of the renormalization of a set of running variables. This 
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is usually done in terms of the relevant parameters of the Hamiltonian. 
The difRculty of the problem at hand is precisely the fact that the effec- 
tive Hamiltonian is not simply related to the microscopic one we started 
from. We have seen in the previous section that the logarithmic correction 
originates from an integration over the collective diffusive density modes. 
The propagator of the diffusive mode, which we have called ladder, is the 
result of a re-summation of impurity scattering to all orders. If we are 
going to obtain a renormalized theory, we must understand how the loga- 
rithmic correction we have obtained in the previous section will manifest 
in the effective field theory for the diffusive modes. This task can actually 
be carried out systematically from the outset by deriving a matrix-field 
non- linear cr- model |28| . This approach, however, makes the connection 
with the physical quantities less direct. By the use of the Ward identities, 
all the renormalizations required by the non-linear cr-model are expressed 
through the Landau parameters, thus showing that the disordered inter- 
acting electron system can be effectively described by a scale-dependent 
Landau Fermi-liquid theory [211 IS E^l EH [73 ■ For this reason, we adopt 
here the approach of perturbation theory combined with the Ward identi- 
ties. 

We begin by recalling the expression of the density and current in linear 
response regime 

where = and = (0, A). iiT^i^ is the standard response function 
for density and current. In particular the conductivity and the compress- 
ibility are given by 

lim , 

uj^Q UJ 

lim Xoo(0,q), 

where K = (1 / d) J^i d assuming spatial isotropy. In a similar way, one 
defines the spin-spin x(w,q) and energy-energy XB(w,q) response func- 
tions, whose static limit (i.e., w ^ first, q ^ after) are the spin 
susceptibility and the specific heat times the temperature. 

The global conservation law and gauge invariance give rise to the fol- 
lowing Ward identities 

qt,Kfj,i, = 0, 

K^,q, = 0, (17) 

that must be obeyed to all orders in perturbation theory. Notice that 
the first equation is nothing but a different way of writing Eq. (Q. For 



a 
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instance, a disordered non-interacting Fermi system has a density-density 
response function given by 

dn Dq^ On iujNq 

9/i Dq'^ — iuj dfi Dq^ — \uj ' 

The first expression is nothing but a statement about the diffusive nature 
of the density fluctuations. The second expression, which separates the 
so-called dynamic part of the response function, is the result of a direct 
diagrammatic calculation. To understand this point, recall that the eval- 
uation of the bubble-like diagram involves an integral over the energy of 
a product of two Green's functions. The dynamic part corresponds to the 
energy range where the two Green's functions have poles on opposite sides 
of the real axis, in the complex plane, and a ladder resummation can be 
performed. A'o appears since we are dealing with a non-interacting system. 
From Eqs. (fT7|l . it follows that 

ifoo(^,0) = 0, (18) 

which immediately gives that the compressibility dn/dn coincides with the 
single-particle density of states, as it is expected for the non-interacting 
system. The single-particle density of states, furthermore, may be generally 
related to the dynamic part of the response function A'^Jq", i.e., 

X„+- = -i /" ^ V G(p, £)G(p + q, e Lo)K+- (p, q, e, lo), (19) 

where is the density vertex in the energy range e(e + w) < 0, which 
in the non-interacting case gives the total dynamic part of the response 
function. The reason for defining the vertex Kq~ in this energy range is due 
to the possibility of inserting the ladder resummation. For positive external 
frequency w, the energy restriction determines the integration range over 
the energy e in Eq. (|19l) . By using the general Ward identity Eq. in 
the dynamic limit, one gets the single-particle density of states 

iV=-- VlmG(p,£ = iO+)= lim lim A:+-(cj, q). (20) 

P 

In the simple non-interacting case, of course, N — Nq. 

We stress that both Eqs. H18I20|I are always valid and we use them in 
the following when dealing with interacting electrons. Here, they will al- 
low us to control the logarithmic corrections which arise from the interplay 
between the diffusive motion of the electrons and their mutual interaction. 
We refer to renormalized quantities with a bar. The diagrammatic skeleton 
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structure of the response function is shown in Fig. is, in the interact- 
ing case, the vertex which, when muhiphed by Kqq~ gives the total dynamic 
part of Kqo^ which includes also terms ending with two advanced (H — h) or 

retarded ( ) Green's functions [37] ■ A direct perturbative evaluation of 

the density response function gives 



where 



Zs = Z- 27VoC'f , (22) 

is the renormalized form of Z^, defined in Eq. and the singular 

corrections have been absorbed in the renormalized quantities D, Z, ^, 
Ts and A^. The first three are obtained by considering the corrections to 
the ladder in the presence of interaction, f ^ is the renormalized scattering 
amplitude in the singlet channel. 

Here, the parameter plays the role of the wave-function renormaliza- 
tion of the effective field theory with the ladder as a propagator, while Z 
renormalizes the frequency, and it will be identified with the specific heat 
renormalization parameter. By the use of Ward identities we now show 
that ^ also coincides with the renormalization of the single-particle density 
of states introduced by the disorder, in the presence of the interaction, i.e., 

N = CA^o. (23) 
Indeed, when Eqs. H18I20|I are used in Eq. (|21|l . we obtain 

No Zs dfi Zs 

Since the compressibility has no singular corrections, Zg coincides with Z^ 
and CAs — whence Eq. H23|l follows. Thus, the single-particle density 
of states becomes scale-dependent, in contrast to the non-interacting case. 

In a similar way, considering Zt ~ Z — 2No('^rt, one obtains for the 
specific heat and for the spin susceptibility the relations 

Cy = ZCy , 
X = ZtNo. (25) 

The above equations relate the renormalization parameters Z, and Zt that 
appear in the perturbative correction to the physical quantities Cy and 
X- Contrary to Zg, Z and Zt are affected by logarithmic singularities, 
which can be derived through the evaluation of the thermodynamic poten- 
tial. Alternatively, one can compute the corrections to Z, and Zt directly 
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(although with a lot more effort) in perturbation theory and find a per- 
fect agreement between the expressions obtained via the thermodynamic 
potential and Eqs. (j^SJl. 

Before going to the next section, to deal with the RG equations, we 
would like to recall the complete expression for the electrical conductivity. 
To do so, we notice that in the presence of Coulomb long-range forces, to 
avoid double counting, one has to subtract the statically screened long- 
range Coulomb Tq from the full singlet scattering amplitude entering the 
ladder resummation for the density-density response function. Hence, f s 
fs — Ea. H22l) for Zg is modified accordingly, and Tq reads 

ro(.,c. = 0) = ^^('^^^^ .,^0^^. (26) 

^ ^ l + Vc{q)2dn/dfi ' 2dn/dn ^ ' 

As a consequence, by using Ea. H24(l . we derive the constraint 

Z - 2NQCfs. (27) 

Then the EEI correction to the conductivity by including also the effect of 
the triplet channel and of the dynamical re-summation |28[ |2] becomes 



5t = -t 



l + 3(l-il±Min(l + 27.; 
27* 



In(Tr), (28) 



where No^'^Tt/Z = 74. The first term in the square brakets is the contri- 
bution due to the singlet channel in the presence of long-range Coulomb 
interaction [see Eq. ||SJ)]. Notice that the corresponding scattering ampli- 
tude does not appear explicitly, due to the constraint (j?7jl . An analogous 
contribution would come from the WL correction, given by Eq. ((TJ. This 
last term, together with the contribution of the interaction in the particle- 
particle channel, is here suppressed to maintain the discussion simple while 
keeping all the general features of the theory. 

Finally, we would like to remark about the various "density of states" 
we have encountered. We have the thermodynamic density of states or 
compressibility dn/dn, the coefficient of the linear-in-T term of the specific 
heat, 7, and the single-particle density of states obtained from the Green's 
function, N. In the non-interacting case all these quantities, together with 
the spin susceptibility, coincide with Nq, but in the presence of interaction 
and disorder they are renormalized by Zs, Z, C, and Zt, respectively. The 
Ward identities are a very powerful tool in controlling the different role 
played by these quantities in the renormalized theory. 
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4.4 One-loop RG equations for the renormalized Fermi 
liquid 

In the previous section we have seen that all the divergences which appear 
in perturbation theory may be absorbed into the renormalization of the 
physical parameters characterizing a Fermi liquid. Therefore, we can use 
the previous results to obtain here the RG flow equations for Z = Cy jCy^Q 
and ZtjZ — x/Zxo- From Eqs. (|15I28|I . introducing the scahng variable 
s = -ln(rT), one gets [HI Ell 132 El EH 

dZ tZ , 

S - i(l + 27*f, (29) 



At t 



ds 2 



l + 3(l-(i±Mln(l + 2,.) 



where the last equation is derived through Eq. I|28|) . considering the bare 
dimension e = d — 2 of the coupling t in units of inverse length. We briefly 
review the main consequences of the above equations. Let us begin the 
discussion in two dimensions, i.e., e = (see Fig. I1U|I . Under the RG, 7* 
grows always and diverges at a finite scale. In the equation for t the sin- 
glet contribution makes it larger (localizing character) , whereas the triplet 
contribution (in the round brackets) does the opposite (antilocalizing). If 
the initial value of 74 is not too large, t initially grows until the growth of 
7t makes the triplet contribution the dominating one. As a result, t has 
a maximum as a function of s. However, due to the strong-coupling run- 
way of 7t, one cannot seriously trust the above equations quantitatively. 
Nevertheless, the physical indication of some type of ferromagnetic insta- 
bility is rather clear due to the diverging spin susceptibility associated with 
7f. Furthermore, the dominating antilocalizing effect of the triplet while t 
remains finite, strongly supports the possibility of a metallic phase at low 
temperature, in contrast with the non-interacting theory based on WL only 
\6'2\ 1341 1571 1381 139| . Indeed, this metallic phase in d = 2 has been recently 
observed (See Refs. in IS?). In addition, due to the divergence of 74 also 
Z goes to the strong coupling regime, leading to an enhancement of the 
specific heat, which is however hardly observable in two dimensions. 

In three dimensions (see Fig. Illll one has a richer scenario depending on 
the initial values of the running couplings. The main difference is that there 
exists a critical line in the t — plane asymptotically given by t^t = e/2, 
where under the RG flow t ^ 0, 7* — > cxd with their product being constant. 
On the weak-disorder side of this line, the system scales to a conductor 
with vanishing t, whereas on the strong-disorder side the system behaves 
qualitatively as in the two-dimensional case discussed above, leaving again 
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the possibility open for a ferromagnetic instability at a finite scale. In 
this latter case, however, the strong-coupling runaway flow requires to go 
beyond the one-loop approximation we have presented here leaving the 
problem of the proper treatment of the metal-insulator transition still open. 
An approximate treatment of the two-loop correction is possible, but its 
discussion is well outside the scope of this paper. We refer the reader to 
Ref. JS]. However, we remark that, whenever there is a reduction of 
the effect of the triplet channel, due to a magnetic coupling altering the 
symmetry, as in the presence of a magnetic field with Zeeman spin-splitting 
spin-flip scattcring[2niEllSSl or spin-orbit scattering one 
obtains a bona-fide MIT with different universality classes with respect to 
the Anderson localization. For a comparison with the experiments see, e.g., 
Ref. [ZH. 

We finally emphasize that, as in the interacting boson case, the pertur- 
bative RG at one loop confirms the conclusions drawn by the implemen- 
tation of the Ward identities. However, in the present case, since we are 
dealing with critical behavior, the singularities are resummed in power laws, 
rather than to be cancelled as required in a stable liquid phase. Mathemat- 
ically this manifests in the lack of additional symmetries which prevents us 
from closing the equations of motion and solving the problem exactly. We 
have achieved, however, via the Ward identities, that the remaining singu- 
larities describe the critical behavior of the specific heat, spin susceptibility, 
and of the coupling t related to the conductance. 
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Table captions 

Tab. ^ In the table we give a summary of the conservation laws 
exploited in order to keep under control the singularities arising in pertur- 
bation theory, for each of the physical systems discussed in the text. In 
particular, we indicate the conservation laws (second column), the practical 
consequence of their implementation (third column), the behavior of the 
system (fourth column), and the equation number where the corresponding 
Ward identity can be found in the article (fifth column). 
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Physical 
System 


Conservation 
Law 


Consequence 


Behavior 


Ward Iden- 
tity 


Interacting 
Fermions 


Global and 
separate Left 
and Right 
number 
density 


Cancellation 
of singu- 
larities in 
correlation 
functions 


Unusual 
non- 
critical 


Eqs. HEl) 


Interacting 
Bosons 


Global 
number 
density and 
t/(l)-gauge 
invariance 


Identification 
and reduc- 
tion of 
running 
variables 


Unusual 
non- 
critical 


Eq. ® 


Disordered 
Interacting 
Electrons 


Global 
number 
density and 
C/(l)-gauge 
invariance 


Identification 
and reduc- 
tion of 
running 
variables 


Critical 


Eqs. piTTjl 



Table 1: C. Di Castro, et al. - Renormalization group... 
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Figure captions 

Fig. n (a): Lowest-order correction to the coupling constant of the (/j*^ 
theory. Each hne in the loop carries a factor [t + q^]~^, and an integral 
over d'^q is implied, giving a contribution which diverges as 1^1"*^^^ for 
t 0, when e > 0, i.e., d < 4. The divergence is logarithmic for e = 
{d — 4). (b): Lowest-order correction to the coupling constant of an 
interacting Fermi system. Each line in the loop carries a factor [e — tii?(|p| — 
Pp)]~^, where vp the Fermi velocity and pp the Fermi momentum, and 
an integral over ded'^p is implied, giving a logarithmic divergence in the 
external momentum or frequency, when d — 1. (c): The four-point vertex of 
an interacting Bose systems generates also two- and three-leg vertices below 
the condensation temperature. The dashed line represents a condensed 
boson, (d): Lowest-order correction to the single-particle propagator of 
an interacting Bose system, associated with three-leg vertices below the 
condensation temperature. Each line in the loop carries a factor [o;^ — 
c^q^]~^, and an integral over dujd'^q is implied, giving a divergence as s"*^, 
where s is a finite external momentum or frequency, for e > 0, i.e., c? < 3. 
The divergence is logarithmic for e = {d ~ 3). 

Fig. [3 Dyson equation for the electron Green's function G (thick solid 
line). The thin solid line represents the bare Green's function Gq, the tick- 
dashed line is the effective interaction D (dressed by the RPA series), and 
the black triangle represents the irreducible scalar vertex Aq. 

Fig. O Left: Various four- and three-leg vertices arising in the broken- 
symmetry phase. Dashed and solid lines represent longitudinal {L) and 
transverse (T) modes, respectively. Right: Equation for the most singular 
part of the two-point vertex vll (filled square). The empty square repre- 
sents the bare vertex vlj^^, whereas the empty and black circles represent 
the three-leg vertices v^rpj, (bare) and vltt (dressed), respectively. The 
external dashed lines, which represent longitudinal fluctuations, are am- 
putated, and are only drawn to indicate the corresponding ingoing and 
outgoing momenta. The solid line represent the propagator of transverse 
fluctuations, Gtt- Notice that the perturbative one-loop RG equation for 
Vll is obtained by taking the bare vertices only. 

Fig. 0] Left: Diagram giving the WL correction to the conductivity, 
Eq. 0. The dashed line represents the average of two impurity insertions, 
introduced later on in the text, Eq. (|10|l . In general, it is necessary to sum 
up diagrams with an inflnite number of impurity lines. The leading term 
gives rise to the ladder diagram explained in the text [cf. Eq. (|12|) , see also 
Fig. [7]. In the case of electrical conductivity, one can show that this leading 
term reproduces the semiclassical result of the Drude-Boltzmann theory. 
The next correction is obtained by considering the scries of maximally 
crossed ladder diagrams shown here. In Rcf. 64 , it has been shown 
that the crossed ladder may be evaluated in terms of the direct ladder by 



34 



exploting the time reversal symmetry. Notice that, because of the crossing 
of all the impurity lines, the impurities visited by the top Green's function 
line are met in reversed order by the bottom one. This is the diagrammatic 
representation of the interference of time-reversed trajectories. Right: The 
WL correction due to the interference by pairs of trajectories, one the time 
reversed of the other, is shown in a more pictorial way. In the figure, the 
two trajectories (solid and dashed lines) differ in the way one goes around 
the closed loop. 

Fig. |S1 Left: After non-trivial cancellations, as explained in Ref. 
the diagrams (a) and (b) are those responsible for the EEI correction to 
the conductivity, Eq. ©. The diagrams shown here contribute to the 
singlet channel. Another set of diagrams, not shown here, may be obtained 
by considering the Hartree-like contributions. This latter set of diagrams 
gives rise to the correction to the conductivity in the triplet channel. In the 
figure, the tick dashed line represents the electron-electron interaction and 
the shaded parts the insertions of the infinite series of ladder diagrams [ see 
also Fig. [7]. Right: A physical picture of the processes taken into account 
by the diagrams on the left. Two trajectories (A and B) that differ by a 
closed loop do not, in general, interfere. The electron-electron interaction, 
however, may lead to interference, since the extra phase gained in the closed 
loop may be canceled by another electron (C) going around the same closed 
trajectory. 

Fig. Left: Self-energy in the Born approximation. The dashed line 
represents the average of two impurity insertions, Eq. When the 

internal Green's function (solid line) is replaced with the dressed Green's 
function one obtains the self-consistent Born approximation. Right: First- 
order correction to the thermodynamic potential (exchange term). Upon 
impurity averaging, the ladder re-summation appears, which corresponds 
to repeated independent scattering events (see Fig. O 

Fig. \7\ The diagrammatic series of repeated impurity scattering. The 
dashed line represents the average over the impurity strength distribution. 
Physically, both the top (electron) and bottom (hole) Green's function lines 
visit the same impurity site. Notice that, since the impurity scattering is 
elastic, the energy is not changed along an electron line. The infinite series 
of ladder diagrams is indicated with a shaded rectangle. 

Fig. IHl Lowest order diagrams for small (Fi) and large (r2 ) momentum 
transfer (scattering angle). 

Fig. 1^1 Skeleton diagram of the response function. The black triangles 
represent the vertex A^, the open squares the scattering amplitudes r^^t, 
and the shaded rectangle the ladder of impurity lines, Eq. H12II . Kqq is the 
static limit of the response function, i.e., the compressibility. 

Fig. El RG fiow corresponding to the equations (j^ . for d = 2 (e = 0), 
with 7 = 7t. 7/(1 -|- 7) is plotted as a function of t/(l -|- t) to transform 
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the interval [0, +00) into the interval [0, 1] on both axes. Observe that 74 
always scales to strong coupling, i.e., 7/(1 + 7) — > 1. 

Fig. ^2 RG flow corresponding to the equations H29|l . for d = 3 (e = 
1). The notations are the same as in Fig. ^| Observe the two possible 
asymptotic behaviors: on the weak-disorder side (small t), 7 scales to a 
finite value, whereas a two-dimensional-like behavior, with 7 — > cx), is found 
on the strong-disorder side (large t) . The separatrix is marked by a thicker 
line. 
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Figure 1: C. Di Castro, et al. - Renormalization group. 
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Figure 2: C. Di Castro, et al. - Renormalization group... 
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Figure 3: C. Di Castro, et al. - Renormalization group. 
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Figure 6: C. Di Castro, et al - Renormalization group... 



42 



e+w 

p+q p'+q 




= L ( q , CO ) 



P P' 

E 

Figure 7: C. Di Castro, et al. - Renormalization group... 
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Figure 8: C. Di Castro, et al. - Renormalization group. 
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Figure 9: C. Di Castro, et al. - Renormalization group... 
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Figure 11: C. Di Castro, et aJ. - Renormalization group... 
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